clear;clc;
% 得到颜色直方图
pictureNum=1000;
channel=3;
colorMatric=[];
for i=0:pictureNum-1
    imgname=['E:\XDU-6\多媒体数据处理\third\corel\corel\',num2str(floor(i/100)),'\',num2str(i),'.jpg'];
    img=imread(imgname);
%     img=img(:,:,channel);
    img=rgb2gray(img);
    [count,bins]=imhist(img);
    colorMatric=[colorMatric;count'];
end

%%
% % 统计准确率
% acc=zeros(10,1);
% trueNum=zeros(10,1);
% for i=0:pictureNum-1
%     [indexs]=pictureSearch(colorMatric,i);
%     pakageIndex=floor(i/100);
%     for j=1:10
%         if floor(indexs(j)/100)==pakageIndex
%             trueNum(pakageIndex+1)=trueNum(pakageIndex+1)+1;
%         end
%     end
% end
% acc=trueNum/1000;

%%
% 实现图像检索
imgIndex=input('请输入检索相似图的图片index\n');
[indexs]=pictureSearch(colorMatric,imgIndex);
% 画图
figure(1);
imgname=['E:\XDU-6\多媒体数据处理\third\corel\corel\',num2str(floor(imgIndex/100)),'\',num2str(imgIndex),'.jpg'];
imshow(imgname);title(['原始图像',imgIndex]);
figure(2);
for i=1:10
    subplot(2,5,i);
    imgname=['E:\XDU-6\多媒体数据处理\third\corel\corel\',num2str(floor(indexs(i)/100)),'\',num2str(indexs(i)),'.jpg'];
    imshow(imgname);title(['图片的index：',num2str(indexs(i))]);
end
